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ABSTRACT 

Information  gleaned  from  applying  the  frozen 
stress  photoelastic  method  to  cracks  emanating  from 
critical  locations  around  the  fin  tips  in  models  of  solid 
rocket  motors  is  reviewed  and  assessed  together  with 
new  experimental  results.  The  studies  are  the  ini¬ 
tial  part  of  a  program  developed  to  contribute  back¬ 
ground  data  for  consideration  in  modifying  current 
motor  grain  design  rationale. 

INTRODUCTION 

Computational  analysis  and  two  dimensional  ten¬ 
sile  tests  on  single  motor  grain  fins  suggest  that 
cracks  in  fin  tips  are  most  likely  to  originate  at  the 
point  of  coalescence  of  the  fin  tip  end  radius  with  a 
small  radius  connecting  the  fin  tip  end  radius  with 
the  side  of  the  fin.  We  call  such  cracks  off-axis  cracks. 
Designers  of  motor  grain  geometry,  however,  assume 
the  critical  crack  to  originate  on  the  fin  axis  as  a 
Class  1  crack.  Frozen  stress  experiments  show  that  • 
the  former  type  of  crack  (i.e,  off-axis)  is  initially 


a  Class  2  crack  under  mixed  mode  load  and  must 
turn  to  eliminate  the  shear  modes  before  becoming 
a  Class  1  crack.  CottereU,  (1995)  in  differentiating 
between  Class  1  and  Class  2  cracks,  noted  that  all 
Class  2  cracks  become  Class  1  cracks  if  allowed"  suf¬ 
ficient  growth.  In  the  present  discussion,  we  shall 
find  that  Class  2  cracks  are  those  sustaining  mixed 
mode  loads  and  will  reorient  their  growth  paths  to 
■eli-minate-the-shear-  modes-  in  order-  to  -attain  Glass-1 
status,  which  exhibit  no  shear  modes  in  our  current 
test  program. 

THE  EXPERIMENTS 

Recent  frozen  stress  experiments  have  been  con¬ 
ducted  by  Smith,  Constantinescu  and  Liu  (2002)  on 
off-axis  cracks  emanating  from  the  coalescence  of  the 
large  fin  tip  radius  Rll  (Fig.  I),  "and  the  small  ra¬ 
dius  R1.3  (Fig.  1)  which  also  coalesces  with  the  side 
of  the  fin.  Starter  cracks  were  initiated  in  two  ways': 

i.  Normal  to  the  fin  surface 

ii.  Parallel  to  the  fin  axis 
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Both  cracks  were  inserted  by  striking  a  sharp  blade 
held  against  the  fin  surface  at  the  above  noted  crit¬ 
ical  point.  The  starter  cracks  emanated  from  the 
tip  of  the  blade  into  the  material.  The  test  models 
were  then  heated  above  critical  temperature  (Ap¬ 
pendix  A)  and  then  internal  pressure  was  applied 
in  order  to  grow  the  crack  to  its  desired  size  after 
which  the  pressure  was  reduced  to  stop  the  crack 
growth  and  the  models  were  slowly  cooled  to  room 
temperature  and  the  load  removed.  Thin  slices  were 
then  removed  mutually  orthogonal  to  each  flaw  bor¬ 
der  and  its  surface  and  analyzed  photoelastically  as 
two  dimensional  models  but  containing  the  three  di¬ 
mensional  effect.  The  photoelastic  data  were  then 
converted  into  stress  intensity  factors  ( SIFS )  using 
a  two  parameter  algorithm  (Smith  and  Kobayashi, 
1993).  The  Mode  I  and  Mode  II  algorithms  for  these 
calculations  are  given  in  Appendix  B. 

For  the  off-axis  cracks  normal  to  the  fin  surface^ 
Fig.  2  shows  two  such  cracks  which  can  be  regarded 
as  the  same  crack  at  different  stages  of  its  growth. 
During  the  first  stage  (Model  1,  Fig.  2a)  the  crack 
turns  all  along  its  border  under  the  action  of  both 
Modes  I  and  II  everywhere  except  at  the  fin  sur¬ 
face  where  only  Mode  I  exists.  Section  SS  suggests 
that  the  starter  crack  was  planar  until  turning  began. 
Model  2}  (Fig.  2b)  shows  another  crack  which  has 
experienced  additional  growth.  After  turning  under 
Modes  I  and  II,  this  crack  exhibits  river  markings 
indicating  Mode  III  as  well  as  I  and  II  until  it  com¬ 
pletes  its  turning  and  becomes  parallel  to  the  fin  axis. 
No  calculations  were  made  near  the  river  markings 
because  small  deformation  theory  may  have  been  ex¬ 
ceeded  there.  At  the  stage  of  growth  shown  in  Model 
2  (Fig.  2b)  only  pure  Mode  I  exists  all  along  the 
crack  border  and  it  has  become  a  Class  1  crack. 

For  off-axis  cracks  parallel  to  the  fin  axis,  a  small 
damage  zone  occurred  when  the  crack  entered  the 
material  and  it  then  emerged  in  a  direction  approxi¬ 
mately  parallel  to  the  fin  axis.  (Section  SS  of  Model 


3,  Fig.  3).  The  crack  shown  in  Fig.  3  had  also 
achieved  Class  1  status. 

Although  ample  evidence  appears  in  Fig.  3  of 
the  presence  of  Mode  II  turning  and  Mode  III  river 
markings,  the  effects  here  produce  what  appears  as 
only  negligible  turning  as  shown  by  Section  SS. 

It  is  this  type  of  crack  which  we  wish  to  make 
comparisons  with  symmetric  cracks  on  the  fin  axis, 
as  shown  on  Fig.  1.  For  this  purpose,  a  series  of 
tests  on  models  containing  a  single  symmetric  crack 
as  located  in  Fig,  1  were  tested  using  the  same  pro¬ 
cedures  as  were  described  above.  To  date  two  such 
tests  have  been  completed  (Models  4  and  5).  One 
such  crack  (Model  4)  is  shown  in  Fig,  4,  and  is  a 
planar.  Class  1  semi-elliptic  crack  with  pure  Mode  I 
all  around  the  crack  border. 

RESULTS 

Test  data  and  results,  including  those  for  Models 
1  and  2  discussed  above  are  presented  in  Table  1. 
Models  3  and  4  are  similar  in  most  respects  except 
that  Model  3  was  an  off-axis  crack  with  blade  parallel 
to  the  fin  axis  while  Model  4  was  a  symmetric  crack 
on  the  fin  axis.  Results  from  a  second,  deeper  sym¬ 
metric  crack  (Model  5)  can  be  compared  to  results 
from  the  average  of  two  off-axis  cracks  with  blades 
parallel  to  the  fin  axis  which  were  in  the  same  test 
specimen  (Model  6)  but  were  separated  by  an  un¬ 
cracked  fin  as  per  Fig.  1.  Finally,  Fig,  5  shows  the 
path  of  the  mid  point  of  an  off-axis  crack  (Model  7) 
with  blade  parallel  to  the  fin  axis  which  grew  to  and 
penetrated  the  outer  boundary  of  the  model.  This 
suggests  that,  once  the  off-axis  cracks  have  cleared  a 
region  of  about  half  the  fin  width  around  the  crack 
tip,  they  will  follow  paths  parallel  to  the  symmet¬ 
ric  cracks  to  the  model  boundary.  It  is  also  noted 
that  the  average  a/c  for  two  symmetric  cracks  which 
penetrated  the  boundary  in  other  tests  was  0.79  as 
compared  with  0.77  for  Model  7. 
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Table  1:  Table  1  -  Data  and  Results 


Loads1 

Crack  size2 
(dimensions  in  mm) 

t 

Notations 

i 

Pmax  =  0.049  MPa 

Psf  ~  0.035  MPa 

Model  1 

off-axis  normal  to  surface 
a  =  8.T1  Aa  =  2.18 
c  =  11.15  Ac  =  3.02 
a/c  =  0.78  a/t  =  0.23 

Fi  =  2.09 
F2  =  0.53 

1-  Pmax  =  internal  pressure  to  grow  crack 
psf  =  stress  freezing  pressure 

Pmax  =  0.103  MP& 

psf  =  0.049  MPa 

Model  2 

off-axis  normal  to  surface 
a  =  12.50  Aa  =  3.40 
c  =  21.10  Ac  =  10.40 
a/c  =  0.59  a/t  —  0.34 

2.19 

0 

2.  a  =  crack  depth; 

Aa  =  crack  growth 

c  =  half  length  of  crack  in 
fin  tip  surface 

Ac  =  half  crack  growth  in 
fin  tip  surface 

3*  Pmax  —  0.103  MPa 

psf—  0.049  MPa 

Model  3 

off-axis  parallel  to  fin  axis 
a  =  7.9  Aa  =  2.80 
c  =  13.35  Ac  =  7.75 
a/c,=  8.59  a/t  =  0.21 

2.12 

0 

3.  Fi  —  Ks/mpsfyfim 
i  =  1,  2,  at  maximum  depth 

\/Q  =  approximation  of  elliptic  integral 
second  kind 

pmax  =0,121  MPa 

Psf  =  0.041  MPa 

Model  4 

symmetric,  on  fin  axis 
a  =  8.13  Aa  =  2.80 

c  =  12.58  Ac  -  6.10 
a/c  =  0.65  a/t  =  0.22 

2.21 

0 

Q  =  l  +  1.464(f) 165  f<l 

pmax  0.129  MPl 

psf  =  0.046  MPa 

Model  5 

symmetric  on  fin  axis 
a  -  14.6  Aa  =  5.1 
c  =  23.0  Ac  =  13.0 
a/c  =  0.64  a/t  =  0.39 

1.71 

0 

All  flaws  were  characterized  as  semi-elliptic  flaws 
of  depth  a  and  length  2c. 

However*  off-axis  cracks  were  neither  perfectly 
semi-elliptic  nor  planar. 

Pmax  =  0.103  MPa 

psf  =  .049  MPa 

Model  6 

off-axis  parallel  to  fin  axis 
a  =  14.75  Aa  =  7.03 
c  =  22.55  Ac  =  13.87 
a/c  =  0.65  a/t  =  0.40 

1.73 

0 

Pmax  ~  0.110  MPa 

psf  —  0.049  MPa 

Model  7 

off-axis  parallel  to  fin  axis 
a  =  37.08  Aa  =  30.9 
c  =  48.0  Ac  =  41.65 
a/c  =  0.77  a/t  =  1.00 

- 
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DISCUSSION 

On  the  basis  of  the  results  to  date,  there  ap¬ 
pears  to  be  a  substantial  similarity  between  the  off- 
axis  and  the  symmetric  cracks,  both  with  respect 
to  normalized  SIFS  (Fi)  and  also  crack  paths,  par¬ 
ticularly  when  the  former  achieves  Class  1  status. 
However,  there  is  some  initial  reorientation  of  the  off- 
axis  cracks  to  eliminate  shear  modes  before  achieving 
Class  1  status,  after  which  they  travel  parallel  to  the 
fin  axis.  Once  this  reorientation  occurs,  it  may  be 
unimportant  to  differentiate  betwen  the  two  types  of 
cracks.  However,  for  the  full  growth  period,  one  must 
recognize  that  the  symmetric  crack  exhibits  Class  1 
status  from  its  inception. 

Experiments  on  deeper  symmetric  cracks  are  con¬ 
tinuing  and  will  be  complemented  with  equal  depth 
cracks  which  run  all  the  way  to  the  ends  of  the  motor 
grain  models.  Final  results  will  be  used  to  assess  cur¬ 
rent  design  rationale  which  is  based  upon  the  latter 
type  of  crack.  These  preliminary  results  however, 
suggest  that  part  through  cracks  break  through  to 
the  outer  circumference  of  the  model  before  running 
the  length  of  the  cylindrical  model  and  might  be  a 
candidate  for  replacing  the  through  crack  used  in 
current  design. 
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APPENDIX  A-  Frozen  Stress 
Photoelastieity 

When  a  transparent  model  is  placed  in  a  circu¬ 
larly  polarized  monochromatic  light  field  and  loaded, 
dark  fringes  will  appear  which  are  proportional  to  the 
applied  load.  These  fringes  are  called  stress  fringes 
or  isochromatics,  and  the  magnitude  of  the  maxi¬ 
mum  in-plane  shear  stress  is  a  constant  along  a  given 
fringe. 

Some  transparent  materials  exhibit  mechanical 
diphase  characteristics  above  a  certain  temperature, 
called  the  critical  temperature  ( Tc ).  The  material, 
while  still  perfectly  elastic  will  exhibit  a  fringe  sen¬ 
sitivity  of  about  twenty  times  the  value  obtained  at 
room  temperature,  and  its  modulus  of  elasticity  will 
be  reduced  to  about  one  six-hundredth  of  its  room 
temperature  value.  By  raising  the  model  tempera¬ 
ture  above  Tc,  loading,  and  then  cooling  slowly  to 
room  temperature,  the  stress  fringes  associated  with 
Tc  will  be  retained  when  the  material  is  returned  to 
room  temperature.  Since  the  material  is  so  much 
more  sensitive  to  fringe  generation  above  Tc  than 
at  room  temperature,  fringe  recovery  at  room  tem¬ 
perature  upon  unloading  is  negligible.  The  model 
may  then  be  sliced  without  disturbing  the  “frozen 
in”  fringe  pattern  and  analyzed  as  a  two-dimensional 
model  but  containing  the  three-dimensional  effects. 
In  the  use  of  the  method  to  make  measurements  near 
crack  tips,  due  to  the  need  to  reduce  loads  above 
critical  temperature  to  preclude  large  local  deforma¬ 
tions,  and  the  use  of  thin  slices,  few  stress  fringes  are 
available  by  standard  procedures.  To  overcome  this 
obstacle,  a  refined  polariscope  is  employed  to  allow 
the  tandem  use  of  the  Post  and  Tardy  methods  to 
increase  the  number  of  fringes  available  locally. 

In  fringe  photographs,  integral  fringes  axe  dark  in 
a  dark  field  and  bright  in  a  bright  field. 
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Appendix  B 

Mode  I  Algorithm 

Beginning  with  the  Grifiith-Irwin  Equations,  we 
may  write,  for  Mode  I,  for  the  homogeneous  case, 


Kx 


(2ttt)  2 


fijW+Vii 


where: 


aij  are  components  of  stress 
Kx  isSIF 

r,  9  axe  measured  from  crack  tip  (Fig.  B-l) 
(rfj  are  non-singular  stress  components 
Then,  along  9  —  7r/2,  after  truncating 


max  -  Kl  .  _o  _  Kap 

nz  (S7tr)i  (87rr)i 


(2) 


where: 

t°  =  /(<r°)  and  is  constant  over  the  data  range 
Kap  =  apparent  SIF 

T^ai  —  maximum  shear  stress  in  nz  plane 


Normalizing  with  respect  to  cr, 


i^p 


,  \/8' 


ray 


_/  ”tt  ~  T7 — rr  +  "  -  (  -  J  (3) 

a(7ra)2  ~  1  ~  ' 

where  (Fig.  B-l)  a  =  crack  length,  and  a  =  remote 
normal  stress 
Kap 


i.e 


r(7ra)t 


vs-  V^isli 


linear. 


From  the  Stress-Optic  Law,  r™ax  =  nf/2i  where, 
n  =  stress  fringe  order, 

/  =  material  fringe  value,  and 
t  ~  specimen  (or  slice)  thickness 
then  from  Eq.  2 


Kap  =  =  ^(8nr)i  (j 

where  Rap  (through  a  measure  of  n )  and  r  become 
the  measured  quantities  from  the  stress  fringe  pat¬ 
tern  at  different  points  in  the  pattern. 


In  the  present  study,  instead  of  normalizing 
K  with  respect  to  d^a)1/2,  we  have  selected 
Psf^Jna/Q  as  the  normalizing  factor  where  y/Q  is 
an  elliptic  integral  of  the  second  kind  approximated 
here,  as  shown  in  Table  I.  An  example  of  the  deter¬ 
mination  of  F\  in  Table  I  from  test  data  is  given  in 
Fig.  B-2, 

Mixed  Mode  Algorithm 

The  mixed  mode  algorithm  was  developed  (see 
Fig.  B-3)  by  requiring  that: 

rH™o|(87tTm)1/2  — (Ki  >K2,rm,6m,TiA=  0 

(4) 

which  leads  to: 

<*> 

By  measuring  9°m  which  is  approximately  in  the  di¬ 
rection  of  the  applied  load,  K2/Ki  can  be  deter¬ 
mined. 

Then  writing  the  stress  optic  law  as: 

^max  _  fn  _  % AP 

nz  (8tt r)l 

where  K*p  is  the  mixed  mode  SIF,  one  may  plot 
fr  vs*  v  r/a  as  before,  locate  a  linear  zone  and 

cr  {ir a)  2 

extrapolate  to  r  =  0  to  obtain  K*.  Knowing,  K*, 
K2/Kx  and  9°m,  values  of  Ki  and  K2  may  be  deter¬ 
mined  since 


K*  =  [{Kxsine°m  +  2K2cOS^)2 

+  (K2sine°J2}i  (6) 

Knowing  K*  and  9°m,Kik  K2  can  be  determined 
from  Eqs.  5  and  6.  Details  are  found  in  Smith  and 
Kobayashi  (1993). 


Fig,  B-l  Near  Tip  Notation  for  Mode  I. 


viauv  .  .  • 

Mixed  mode  region  Mixed  mode  region 

magnification  factor:  2.75 

Fig.  2a  Model  1  Off-Axis  inclined  Crack  Showing  Starter  Crack  and  Final  Crack  Front 


C  -  crack  front 

D  -  camera  view  of  the  photograph  magnification  factor:  2.00 
Fig.  2b  Model  2  Off-Axis  Inclined  Crack  Showing  Starter  Cracks  and  Final  Mode  I  Crack  Front 
Figure  2:  Typical  Off-Axis  Inclined  Cracks  (Blade  Held  Normal  to  Fin  Surface) 
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Off- Axis  Crack  Path  (Blade  Held  Parallel  to  Fin  Axis).  Model  7 


